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I. INTRODUCTION 

Quantum electrodynamics (QED) is liable to correc- 
tions if the electromagnetic environment of the system 
under consideration is different from free space. For ex- 
ample, the Lamb shift in an atom changes if the atom 
is located not in free space but near a reflecting sur- 
face 0, . This and similar boundary-dependent effects 
are the subject of cavity QED Under most circum- 
stances cavity QED effects are non-relativistic in nature 
and hence the techniques employed in the theory of cav- 
ity QED are chiefly non-relativistic. They rely mostly 
on a comparatively simple mode expansion of the elec- 
tromagnetic field and on first-quantized theory for the 
remaining part of the system under investigation (cf. 
e.g. jH). However, there are a few examples of systems 
that are not inherently non-relativistic, the simplest be- 
ing a free electron. By virtue of being free it lacks an 
in-built energy scale that could limit its virtual excita- 
tions and thus its interactions with the electromagnetic 
field to non-relativistic energies. Other effects that re- 
quire a fully second-quantized theory are, e.g., radiative 
corrections to the Casimir force between reflecting planes 
m and the Scharnhorst effect of faster-than-c light prop- 
agation inbetween and perpendicular to parallel plates 
[g. To date all such field-theoretical calculations have 
been done for cavity walls that are perfectly reflecting. 
While it is obvious that no real material can ever really 
be perfectly reflecting, the model of perfect reflectivity 
seems to capture all the essential physics of the bound- 
ary without going wrong by anything other than a mi- 
nor numerical prefactor. Its great attraction is that it 
is comparatively simple; for example, the photon propa- 
gator between two parallel perfectly reflecting plates can 
be written as a sum of the photon propagator in free 
space and a small boundary-dependent correction, and 
loop calculations using it are manageable, though not 
trivial because of the loss of translation invariance per- 
pendicular to the plates 0. 

However, we recently discovered that the assumption 
of perfect reflectivity for the cavity walls is in fact not 
justified for systems, such as a free electron, that admit 
low-frequency excitations, i.e. whose excitation spectrum 
has, unlike an atom's, no natural IR cutoff 0. This is be- 
cause, crudely speaking, the electron's virtual excitations 



of arbitrarily long wavelengths interact with evanescent 
electromagnetic field modes originating inside a cavity 
wall, and the chief defect of the perfect-reflector model 
is that it ignores all such evanescent modes. To show 
this we have modelled an imperfect reflector by a non- 
dispersive dielectric and have calculated the self-energy 
of an electron in front of a dielectric half-space. We have 
found that taking the limit of perfect reflectivity in the 
result disagrees with the corresponding calculation that 
assumes a perfectly reflecting wall from the outset. The 
two results differ by a factor of 2 in one direction and 
even by sign in the other. While the effect as such can 
already be seen in a non-relativistic calculation, we felt 
that there was a need for a proper second-quantized cal- 
culation, mainly for three reasons: (i) the non-relativistic 
calculation yields different results for the two models but 
gives no clue as to the origin of this discrepancy; (ii) 
there is nothing that a priori restricts the electron's mo- 
tion to non-relativistic energy scales and, in fact, the in- 
teraction energy is being integrated up to infinity, which 
could potentially lead to errors that pass by unnoticed in 
a purely non-relativistic calculation; and (Hi) there are 
other effects, as e.g. the boundary-dependent g — 2 cor- 
rection to the electron's anomalous magnetic moment, 
which might be affected and whose calculation requires 
field-theoretical methods. Furthermore, in the face of a 
major discrepancy it seems wise to check all possibilities. 

In this paper we establish the major building blocks of 
a full QED theory near imperfect reflectors. We concen- 
trate on a non-dispersive and non-absorbing dielectric 
as a good model for an imperfectly reflecting material. 
Thus the medium is characterized by a single parameter, 
its refractive index n, which is real and the same for all 
frequencies. For technical simplicity we restrict ourselves 
to a single reflecting surface, i.e. we consider a dielectric 
half-space, which we take to occupy the region z > 0, 
while the region z < is vacuum (cf. Fig. QJ. For this 
set-up the dielectric function is a single step function 

e(r) = \ + {n'^ ^ I) e{z) , 

which makes the solution of Maxwell's equations com- 
paratively simple. For piecewise constant dielectric func- 
tions like this, it is advantageous gauge to work in the 
generalized Coloumb gauge 

V • (eA) = , (1) 
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FIG. 1: A sketch of the simplest system in cavity QED, con- 
sisting of a quantum object, e.g. an atom or an electron, that 
is located a distance a away from a single reflecting wall. If 
the wall is perfectly reflecting then the electromagnetic field 
is excluded from the region 2: > 0. A more realistic but still 
simple model is to describe the wall by a constant, frequency- 
independent refractive index n. 



which we shall do in this paper. For a general coordi- 
nate dependent dielectric function e(r) the generalized 
Coulomb gauge may be a very awkward choice, but for 
a piecewise constant e(r) this gauge is so convenient be- 
cause anywhere except right on the boundary (or bound- 
aries) of the dielectric {z — in our case), it is equiv- 
alent to the Coulomb gauge V • A = 0. Thus one can 
work almost as if in the Coulomb gauge and just needs 
to make sure that the physical fields satisfy the appro- 
priate matching conditions at the boundary, i.e. that 



Ell continuous , Dj_ continuous 



(2) 



Since our model material is just a dielectric and has 
a magnetic permeability /i = 1, the magnetic fields 
strengths B and H are also continuous everywhere. We 
note that quantum mechanics in the generalized Coulomb 
gauge is related to that in the true Coulomb gauge by 
canonical transformation the two differ by surfaces 
charges leading to an electrostatic image potential in the 
Hamiltonian. 

If one wanted to work in a gauge that resembles the 
radiation gauge, one ought to choose the gauge 

In this gauge Maxwell's equations for the scalar potential 
$ and the vector potential A read 



V(V.A)-V^A + .^-eV^^ 



It is important that they separate and that for piecewise 
constant e(r) they differ from the standard wave equa- 
tions for $ and A only by surface terms. However, in 
the present paper we shall not pursue any calculations in 
this gauge but work in the generalized Coulomb gauge 

In Section m we calculate the full photon propagator 
in the presence of a non-dispersive dielectric half-space 
starting from the normal modes of the radiation field 
which we briefly discuss in Appendix^ In Section lllll we 
use this photon propagator to determine the self-energy 
of a free electron located outside and a distance a away 
from the dielectric. In Section Hvl we do a careful asymp- 
totic analysis of the expression for the self-energy for non- 
relativistic mean energies. In Section ITvnl we compare 
the results for the electron's radiative self-energy in front 
of an imperfectly or of a perfectly reflecting surface and 
discuss the reasons for the disagreement between the cal- 
culation for a non-dispersive dielectric and for a "perfect 
reflector". Section Ivl summarizes our final results. 

Throughout this paper we set c — \ ^ h and use 
Heaviside-Lorentz units for electromagnetic quantities, 
eo = 1 = Mo- Thus the fine structure constant is 
eV47r w 1/137. 



II. CALCULATION OF THE PHOTON 
PROPAGATOR 



A. Wightman function 

The Green's functions in free quantum electrodynam- 
ics are vacuum expectation values of products of field 
operators. Let us first consider the Wightman functions 



op 



D^^^'-ix^x') = -i < 0\A^'{x)A''{x')\Q > 



(3) 



,V • (6A) 



Inserting the normal modes (|A2(I and l|A5(l into Q and 
taking the vacuum expectations values of the bilinear 
products of photon annihilation and creation operators, 
we obtain 

Df^-^-ix^x') 

x[^D;ix,x') + ^'D;{x,x')]. (4) 

with 
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+e{-z)e{z')- U^k.-<-')TR + e'k.-+^kiz')J.RJ^R*^ 

+e{z)e{-z')- \e«-'K-')Tf + e-^'='-»'=:-')i?fl.rfl.*l \ ^ (5) 



and 



<7 (7 



+0{-z)e{~z')[e'''^^'--'-''> + e'''^^'-+''''>R^{k)* + e-'''^^''+'-''>R^{k) + e-'''^^^^ 

+9(-z)6(z')-\e^'''^'-'''"''^T^* + e(-'k^'-^ktz')J.L* j^Li 
n 

+ 0{z)e{^z')-\e^'''''-'''''">T^ + e('ki.+rk^z')j,Lj^L*l 1 

n J 



(6) 



The notations for the wave vectors is as defined in Eq. IIA4p . and in addition we have introduced the new variables 
k = (fc'^,k||) and x = (a;'^,X||) in 2+1 dimensional Minkowski space. The sum of '"D^{x,x') and ^D^ {x,x') can 
be simplified by taking into account that the Fresnel coefficients are real functions of the wave vectors and by using 
various relations (|A7|) and ljA8|l between them and their products. We obtain two equivalent expressions: 



D,{x,x') + ''D^{x,x') = I 



(Pkii 1 



-i{x — x')k , 



27r)3 2uj 

+e{-z)e{-z') 



{z)e{z')- / dk 



+e{-z)0{z')- / dfc^e'*=^^-''='^'rj^ + 6l(z)0(-z')- / dk^e-'^'''+'^''-''T^' 



(7) 



and 



D^{x,x') + "D„{x,x') = I 



d^fcii 1 



(27r)3 2w 
+9{~z)e{~z') 



— OQ 

oo 



+ / dk,e'''^^^+''^ Rj; 



1 r 1 r 

+e{-z)e{z')- dkie'^'''-^'i'''T^ + e{z)e{-z')- dk^e-'^^'^'+^'^Tf \, (8) 



with r = ((n^ — f)fcy)^/^. The difference between these 
two expressions is how the contributions from evanescent 
waves are included. In the first expression IJJ) they ap- 
pear in the second line as a separate integral over kf from 
— r to 0. In the second expression ((SJ they are included in 
the integration along the path C in the complex kz plane 
shown in Fig.|21 it runs along the real axis from — oo to 0, 
then down the negative imaginary axis from to — zF/ti 
to the left of the square root cut, back up to the origin to 
the right of the cut, and then along the real axis from 



to +00. The cut is due to kj = (n^fc^ + {n^ - l)fc|)i/2/n 
and extends from k^ = +iT /n to k^ — ^iV jn. The part 
of C that runs left and right of the cut is identical to the 
integral over fcf in the second line of (0, i.e. it gives the 
contribution of the of the evanescent waves. This works 
because: 

\k^ = -iK.k^ = -K \k;, = -iK.,k'l=K 

= ^nW* . (9) 

'^z k^ = -iK,k'^=-K 
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the Wightman functions Q , and (jSJ according to 



Rek 



= 9{x^ -x'^) D'"'-{x,x') 
D''^'-{x',x) 



■1 ^ 2 ,,1/2, 

-ik|| (n -1) /n 
-ik„ 

<K +K ) 



Thus D^^'^'^{x, x') can be written in the same way as the 
Wightman function in Q), 



FIG. 2: The integration path C in the complex plane. 

Note that the w = singularity at = ±ifc|| does not 
come into play if one chooses the cut of w = (fc^ + fcy )^/^ 
along —ico ... — ifcj| and i/cji . . . ioo. 

B. Feynman propagator 

For the calculation of radiative corrections we need the 
Feynman propagator, which can be reconstructed from 



D'^'^'^ix^x') = ^e;^(9,)er(5.0 g^" Dl{x,x') . 



The sum over the polarizations is gauge dependent. In 
Coulomb gauge the sum runs over a = TE, TM only. 
In covariant, i.e. Feynman gauge the two unphysical po- 
larizations G, C have to be included. Proceeding from 
the simplified expressions (0 and ((HJ for the Wightman 
functions, we obtain for the polarization component a of 
the Feynman propagator 



D%{x,x') = - 



(Pk 



-e{-z)0{-z') 
-e{-z)0{z')- 



e{z)e{z') 



dkt 



^21,2 jL2 l-d2 I ,v 

t rtQ rip n.^ ~p 



1 



k^ + ie 



dkt 6^'=-^-*'=^^' T? 



— OO 

oo 



1 f 

^e{z)e{-z')- / dki e-'^'''+'^'' T^* 



kl - {k^ + fc^2)/„2 ^ 
In an alternative formulation one can replace the second line of by 



i-z)ei-z') 



dk, e'fe^(^-^') 



dk, e'^-^'-+'"^ + 

) 

I 



fc2 



(10) 



(11) 



Note that the property e'^{dx)e'^ {dx') — e'^* {dx)e'^{dx') This is easily shown by a change of variables fc —k. For 

real k^ this is straightforward; for imaginary k^ — —in 
(with K > 0) one needs the relation 



follows directly from the definition IjAip of the polariza- 
tion vectors. Note also that fco in Eq- l|10|l is a free in- 
tegration variable and is not fixed to lj as it was in the 
case of the Wightman functions {7)) and ||HJ). 



T (A^_2 ) ^ (^2 ^/v, 



The propagator differs from its equivalent in free space 
by not being translation invariant along z. However, it is for which one has to take care to stay on the same sheet 
still symmetric in its arguments, D'^ix^x') = D^{x',x). in the complex plane. 
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The wave equations that the Feynman propagator sat- Checking the wave equation for example for z < leads 
isfies are to 

{dl^ - ^^)D%{x,x') = S^'^Hx -x'), for z < , 
{r?dl^~~ /\^)Dl{x,x') = 5^'^\x~x'), for z > . 



(9^ - ^,)D%{x,x')^e{-z)9{-z')5i^\x - x') 



The second line should vanish, which it does — because 
the integration path can be closed in the lower A:^ plane 
as the Fresnel coefficients do not have poles there. 

A direct check of the Feynman propagator is the re- 
construction of the Wightman functions Q and (jHJ by 
performing the feg integration. For xq > the integration 
path can be closed in the lower half plane 



- exp(-i|k|a^O) e{xa) , (12) 



k^ + k'^ + ie 2|k| 



and we thus reproduce and l|Hl for xq — x'q > 0. 

Next we consider two special cases. The simplest sit- 
uation is the limit rt — > 1, when there is no medium. 
Correspondingly all reflection coefficients R tend to zero, 
the transmission coefficients T approach one, and F — > 0, 
so that complex values do not appear in the integration 
over kz- For all polarizations a the propagator compo- 
nents equal the massless scalar propagator 



-i{x — x )k 



1 



fc2 



(13) 



Thus, according to HA9|I our Feynman propagator l(Tn|) 
coincides in the limit n — > 1 with the standard free 



space propagators: either with the covariant propagator 
in Feynman gauge or with the propagator in Coulomb 
gauge, depending on which modes have been included in 
the sum over polarizations. 

In the limit n ^ oo only the left incident modes sur- 
vive, so that, according to (|A10|I . the Wightman func- 
tions |0Jand its components (U)) simplify greatly. 



{x,x') = i 



d^k\\ 
(2^ 



'^e-''^--^'y^6{-z)e{-z') 
2lo 



^ UkAz-z') _^ e*'=^(^+^')i?^(fc)l (14) 



with 



A corresponding representation follows for the Feynman 
propagator. 

Finally we need to establish the electrostatic Green 
function, which corresponds to the Coulomb interaction. 
One way to start is from the retarded Green function for 
the Coulomb mode C. The retarded propagator differs 
from the Feynman propagator (|10|l by the ie prescription 
in the denominator, such that ie — > ie/co- Since we are 
looking for a static Green's function we need to calculate 



Gc{x,x) ^ / i?c*(^:^') 



(15) 



(2^3' 



pi(x||-X||)-k|| 



k2 



— <^ e{-z) 9{-z') / dkz f,ik.{z-z ) ^ ^ik^z+z ) 



,M 1 - 



+ e{z) e{z') j dk, U^^'^^-^') + e-'k^z+z')^^ — \\ 

/OO r\ nOQ r\ 
dkz e^'^^("-"')- ^ + 9{z) e{-z') / dkz e''=^("-"')- ^ 
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In our conventions A^{x) = $*^(a;) for z < and A^{x) = n^'^(x) for z > 0, so that we obtain for the electrostatic 
Green's function 



G^{x,x') = -i(0|$°(x)$"(a;')|0) 
(2^ 



+ e{z) e{z') \ ( e''=-(^-^') + e-^'=-(^+^')^^J^^ 

\ + 1 y 

+ e{-z) e{z') + 

1 + 



(16) 



1 + 



I 

This result agrees with the classical electrostatic Green's momentum p. The particle eigenspinors Mq''(p) are so- 
function, which can be derived easily, so that we have an lutions of the Dirac equation, 
additional check. 



III. THE SELF-ENERGY OF THE ELECTRON 

The energy shift of the electron can be determined by 
considering the electron propagator and its radiative cor- 
rections due to the coupling to the photon field We 
start be considering the free electron propagator 

5S(a:,y) = (0|TVa(x)V^^(y)|0). (17) 
Using canonical quantization one can write the spinor 



+ m 
2m 



( <t>^'^ 
P 



(19) 



with 4>'^^^ two orthogonal and normalized two-spinors. 
Thus the normalization of u'*' (p) is 



7(0 



(p)^.«(p)=<5. 



The antiparticle eigenspinors u(*-'(p) are similar to H19|l . 
t (i) iBt-ip r] except with upper and lower components interchanged 

~^'^i,p^a (p)^ ^ 1 (18) and normalized to -1. Inserting the canonical expansion 

(|18|l and its conjugate into the expression for the propa- 
where 6,,p annihilates an electron of hclicity i and mo- gator ^ and Fourier transforming to go from the time 
mentum p and creates a positron of helicity i and variable {xq — j/o) to the energy E one finds 

I 



S("'rx V n - - / V "''Hp)^^^'Hp) .p.(x-y) , [ ^jr^ V- ^(')(P)^^'''(P) ^^.p.(x^y) ,30) 

i„,(x,y,i.)- y (2,)3s(p)^Z.^s(p)„^_,,^ +7 (27r)3i?(p)^2^^i?(p)+i?-z.^ • ^'"^ 



This shows that the propagator has positive energy poles at the particle energies £'(p) and negative energy poles at 
the antiparticle energies — i?(p). Radiative corrections represent a perturbation and cause shifts in the particle and 
antiparticles eigenfunctions and in their energies. For small perturbations these shifts are small and expressions can 
be linearized in them. Linearizing the change of the propagator SS = S — S^'^\ one obtains a term that is linear in 
the energy shift and has a double pole at the particle or antiparticle energy, 



\ u(')(p)ll(')(p) zjW(p)t;W(p) _,p.(^_^) 



0/5^ ^ - / {2^fE[^)^^^;y^E{v)-E^^,f [£;(p) + i?-.6f ^ 

(other terms without double poles) . (21) 

I 



Further terms that are linear in the shifts of the particle with double poles since the eigenfunctions appear only in 
and antiparticle eigenfunctions do not give rise to terms 
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the numerator of Eq. H20() . 

The energy shift SE{p) can now be determined by com- 
paring the expression H21|l to the change of the electron 
propagator as determined from standard Feynman per- 
turbation theory. At one- loop level, i.e. to order e^, the 
radiative correction to the propagator is 

,55i°j(x,y,i?)= J d'z I dV5i°)(x,z,i?)I],«(z,z',i?) 

X sfiiz',y,E), (22) 

where Sq^(x, x',i?) is obtained from the standard elec- 
tron self-energy 



by Fourier transformation from the time into the energy 
domain. The electron propagator S'^^\x — x') is the same 
as in free space and thus translation invariant in all 4 di- 
rections, but the photon propagator Z?^^(x, x') is affected 
by the presence of the dielectric medium and therefore 
not translation invariant in the 2:3 direction. Substitut- 
ing the representation (|2l)l) in terms of eigenfunctions into 
Eq. one obtains 



J5S(x,y,i?) 



d^z / d^z' 



dh 



dh 



(2^)37 (2^)3 iJ(p) i5(p') 



„ip-(x-z)+ip'- (z'-y) ^a"* (P)"7 '' (P) 

E{p) - E-ie 



X i;^k(z,z',£;) 



E{p') - E-ie 



(24) 



Further terms all contain antiparticle operators and at 
least one negative energy pole. Since we are interested in 
the energy shift of a particle rather than an antiparticle, 
we need to focus only on terms with two particle poles. 
For p = p' Eq. 124(1 has the same double pole as Eq. 
(|21|l . and thus a simple comparison of the coefficients of 
those double-pole terms should yield an expression for 
the energy shift SE{p) in (|21|l . However, a mathemati- 
cally clean comparison is possible only if one introduces 
a quantization volume with periodic boundary con- 
ditions so as to discretize the momentum p. Then all 
integrals over momenta turn into sums according to the 
prescription 



d^p 
(2^ 



-y 

L3 ^ 



Only the term with p = p' in the double sum over mo- 
menta in Eq. (|24|l gives rise to a double pole in the 
energy, and comparison with the double-pole term in Eq. 
((21|l therefore yields 



6Eip) = j^jd^zjd- 



''z'^e-^P<^-^'^ui;\p) 



Eip) 



X S^«(z,z',i?)4*)(p) . (25) 



It is advantageous to work with the Fourier representa- 
tion of E. Because of the lack of translation invariance 
in the X3 direction, the Fourier transform of the self en- 
ergy with respect to x — x' has a residual dependence on 

X3 + Xg, 

^af3[X,X ) - J 



s^^(q^a.3 + 4) ■ (26) 



In a finite quantization volume the integral over d'^q 
again turns into a sum, and we can re-write Eq. 125(1 
as 



-p)-(z-z') — (i) 



X E^,(g,^3 + 4)4'HP) ■ (27) 

Since we want to work out the energy shift of a particle as 
a function of its distance from the dielectric, we need to 
form localized wave packets in the X3 direction — so that 
the concept of a certain distance between the electron 
and the surface of the dielectric at all makes sense. If the 
centre of the packet is at 0:3 = —a (see Fig. then we 
can approximate Z3 -I- Zg ~ — 2a and carry out the z and 
z' integrations in 1(271) . The result simplifies to 



SE{p) 



Eip) 



7(0 



(p)S^4p,-2a)M«(p) . (28) 



Note that, while S(q, 2:3 -I- ^3) in Eq. 1(26(1 is in general off 
shell, it is on the mass shell in Eqs. ((27(1 and ((281) because 
q and p are on shell. 

Further we need to remark that in Coulomb gauge the 
energy shift is not wholly due to the radiative self energy 
((23(1 : we have to add to 1(28(1 the electrostatic energy 



(^i^Couiomb = Y J J ^' "^^^n i^ix) ■■ 

xG^{x,x') ■.i^{x')-f°ij(x') :) , (29) 

where is the part of the electrostatic Green's function 
((16(1 that depends on the presence of the dielectric. The 
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electrostatic shift is easy to evaluate, which we shall do 
in Section HYH 

We are interested in the self-energy corrections for an 
electron located well outside the dielectric. Because of 
the electron's localization its direct interaction with the 
dielectric medium is completely negligible, i.e. there is 
no wave-function overlap between the electron and the 
microscopic constituents of the dielectric. That is why, 
for X3 < and 2:3 < we can work with the standard 
free electron propagator, 

S'-'>'^{x - x') = {ij^'d,. + m)D^^ix - x') , (30) 



D'^{x-x') = - 



(27r)4 



1 



The impact of the dielectric medium onto the self-energy 
of the electron is consequently just due to the electro- 
magnetic interaction, i.e. due to the fact that the photon 
propagator H10|l depends on the presence and electromag- 
netic properties of the medium. Since we are interested 
only in the energy shift due to the presence of the dielec- 
tric, we split the photon propagator into the free photon 
propagator and a medium-dependent part, 

D'^ix, x') ^ DIqj {x, x') + D{x, x') , 

and take only the medium-dependent part D(x,x') for 
calculating the self-energy and the energy shift 
This also means that we do not have to deal with reg- 
ularization and renormalization; these have been done 
in the free part of the photon field, and we work with 
already renormalized quantities. All medium-dependent 
corrections will then automatically be finite. 

From ((Tn|l we see that for z < and z' < 0, i.e. outside 
the dielectric, the medium-dependent part of the photon 
propagator is 

D^''\x,x%.y<o = -Y.e^,{dx)e:*{dx,)g'^'^Dl{x,x') 



where 111,71 = 1,2. For evanescent waves one has fc^ = 
—IK, k1 — —k'^, and = fcy — k'^. 

Inserting the electron propagator and the photon 
propagator (|31|) into the expression for the self energy 
(|23|l . we have to multiply several 7 matrices. Using 



^A'^^T,- = S"^^''7„ - Ze"^^''757a, £0123 = +1, 



we encounter the following 7 valued invariants 



(33) 



T<7 a^Xv 

^25 — ^ 

ja _ afj.Xv 



with q being the Fourier variable in the electron propa- 
gator H30() and g'^^ having the k dependence as in H32(l . 
In terms of those the distance-dependent part of the self 
energy is 



E = ^ = [E? + EJs + E- + EJ5] (34) 



with 



Dl{x,x') 



(31) 

The derivatives of the polarization vectors IjAip act on a 
plane wave and its reflection. We obtain 



with 



grsik) 
5'tm(^) 



h h 



h h 



-k^ 



9 



9 

.00 



3m 



h h h"^ 










k2 








k2 ' 


9'' = 


k2 J 


h h 




kt 


k2 ' 


9'' = 





(32) 



E^(x,2;') = -ie' 



d^k 



(2^)4 y (2^)4 



dk-. 



-iq{x—x') 



-inq,k) — i?^ 



(35) 



In the same way as for the total self energy in Eq. H26(l 
we perform a Fourier transformation for the components 
in the sum H34() . which again retain a dependence on 
z + z' due to the broken translation invariance in the z 
direction. 



Eg (a;, x') — 



{2n) 



E^(p,z + z') 



Making the variable replacements q = p — k, q^ = Pz in 
Eq. H35|) we obtain 
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Eg (p, z + z ) — —ie^ 



(2^ 



dk. 



{p-kf 



Pi 



1,2 



(36) 



This expression looks much simpler than it is to evaluate. The loss of translation invariance perpendicular to the 
surface of dielectric is one source of complications, and the interference of incident and reflected waves is another. In 
order to evaluate the self energy components (|36|l we need the explicit expressions for the invariants Ig , which depend 
on the mode a. For the two physical modes we find 



tTE _ 
-'2 — 



7|| • (P|| - k||) - 7|| • k| 



k|| ■ (Pit -kjl) 



+ 7|l • (p|| -k||)+73p2 -"foipo-ko), 



(37) 



tTM 



k2 

(7n ■k||) fki. 



h h 



(Pii 



fc2k2 



7|| ■ (Pii - kii) + 73p^ - 7o(po - ko) 



(38) 



IV. ASYMPTOTIC ANALYSIS OF THE 
SELF-ENERGY 

A. General approach and approximations 

Our aim is to determine the energy shift of an electron 
that is localized in z-direction. The shift will depend 
on the distance a of the electron from the surface of the 
dielectric, and without localization the notion of this dis- 
tance would not make sense. Physically the localization 
could be realized be sending a tightly focussed beam par- 
allel to the surface or by confining the electron by means 
of magnetic and/or electric fields. So, we will in fact 
not be working directly with momentum eigenstates H19|) 
but use them to form with wave packets that peak at 
z = a and whose average momentum in z direction is Pz ■ 
The wave packet may move as a whole, which is why we 
have not approximated z — z' in Eq. (|27|l : but we shall 
assume the electromagnetic field to be the same across 
the packet, which corresponds to the dipole approxima- 
tion in atomic physics, and which is why we have set 
Z3 + z'^~ -2a in Eq. if^ . 

The extent of the wave packet must be small compared 
with the distance a from the surface, but otherwise the 
details of the wave packet are not relevant. This im- 
plies am ^ 1, i.e. that the distance a must be very 
much larger than the Compton wavelength Xc = f^l {mc). 
Thus, we shall aim for an expansion in l/{apo)- 

In order to proceed with the calculation of the self en- 
ergy H36|l . we want to perform a Wick rotation fco — ik^. 
By design, the poles of the photon propagator he in the 
right position for this. However, the poles of the term 
that originates from the electron propagator may inter- 
fere; they lie at fco = Po ± [^/(Pll ^ ^11 )^ + + ~ 



There is no problem if they come to lie in the 2nd and 4th 
quadrant of the complex ko plane, but for 2p|[ • ky > ky, 
one of the poles lies in the 1st rather than the 2nd quad- 
rant, if we take p to be on-shell, and a pole in the 1st 
quadrant interferes with the Wick rotation. There are 
several ways of dealing with this problem. One could 
work with a strongly deformed integration path and then 
carry along the separate contribution from the pole; or 
one could go off-shell to move the pole out of the 1st 
quadrant and then do an analytic continuation to a re- 
sult for on-shell p; or one could avoid the problem alto- 
gether by approximating p|| « in the denominator of 
(|36|l . We have decided on the last approach because it 
is straightforward and we are not interested in ultrarela- 
tivistic motion of the electron. 

Thus we set py = in the denominator of H36I) but 
leave untouched elsewhere, i.e. retain it in Ig l|S7|l and 
(|38|l . Then we can perform the Wick rotation fco = ik^ 
without problems and obtain 

/cPkw 1 
(^L*^fc| + 2.pofc4 + k| 

lj^'^hi + \.i + kf'- 

The integration over the three dimensional (Euclidean) 
space (ky , ki) can be carried out in spherical polar coor- 
dinates by defining /c4 = pcosi?, fcj. = p smd cos Lp,ky = 
psini? sin We find 

Y.'^{p,Z + z') = -—r dp d{cOSl^)— -J- 

(27r)3 7o ' ' p^ + 2ipopcos^ 

-jjk.-^R-^-j^ d^n. (39) 
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The only ip dependence is in the invariants /g ; carrying 
out the integration and using the fact that p is on-sheh, 
one gets 



with 



1 

2^ 



dip (ll^ + /J^) = 7oipcosi? - 7„ • p„ , 



27r 





27r 



— / dip {I™ + 1^') = -7o*P cos I? + 7|| • Pll 



k2 



(^27oipcosi9 - 273p^ - 7|| • p| 



rdpi™=Q, rd^i^^o. 

Jo Jo 

The next step in the evaluation of Eq. (|39|l is to carry 
out the integration over by means of contour integra- 
tion. The remaining two-dimensional integral over p and 
cos I? can then be calculated asymptotically for \z + z'\ 
very much larger than the Compton wavelength. Since 
the technical details differ between the TE and TM po- 
larizations, we consider their contributions one after the 
other. While the calculation to follow is perfectly general 
for all values of z and z', provided |z-|-z'|po 3> 1, we now 
simplify the notation and set z ~\- z' = —2a, as this is 
the value at which we need to evaluate the self energy in 
Eq. for the radiative shift. 



B. TE contributions to the self energy 

For the TE polarization the integrand of (|39|l has only 
one pole in the lower kz plane, and that is at kz — —ip = 
— i(k| + (cf. Fig.EJ. The kz integration can thus 

easily be carried out by deforming the contour C and 
evaluating the residue at —ip. We emphasize that when 
evaluating the reflection coefficient R^e^ Eq. (|A3|) . at 
this point, one must take great care that the branch cut 
of the square root in kf is indeed taken to run as shown 
in Fig. 121 Renaming cos?9 — t, we can write the result of 
the contour integration as 



/ dpi dt 



loipt - 7|| • P|| 
p + 2ipot 



27oPoi^ - 7|| • P|| 



RrEi^) - 



1 - ^(n2 - 1)<2 + 1 
1 + - l)t^ 



(40) 



Next we scale the integration variable p — 2pot(^. In 
terms of the new variable ^ the integral reads 



e 

4^ 



d( 



dt t, — 



■ 7|| • Pll 







^2 + 1 

X R^fsit) e-^f""*? . (41) 
This integral can be evaluated asymptotically for large 
values of poa. A standard method of obtaining an 
asymptotic expansion for integrals with an exponentially 
damped integrand is repeated integration by parts. How- 
ever, in two-dimensional integrals like the one above, this 
method generally fails because integration by parts in 
one variable generates inverse powers of the other vari- 
able and the resulting integral diverges at the lower limit. 
For a general discussion of this problem and its remedy, 
we refer the reader to Ref. ^Jl- Here we observe that 
R!f^{t) and hence the integrand of 1)41(1 actually behave 
as 0{t^) for t ^ 0. Thus we can integrate by parts in 
the ^ integral twice without jeopardizing the convergence 
of the t integral. In this way we find to leading order in 
l/(Poa) 



1 2jopot'^ - 7|| ■ 



47r2 (4poa)2 Jo 



dt 



RTE{t) 



(42) 



The t integral in this expression is elementary. 



C. TM contributions to the self energy 



The TM polarization is more difhcult to deal with, 
since the invariants ij^^ introduce a factor 1 /k2 into the 
integrand of (|39() . which leads to an additional pole in 
the lower kz plane at kz ~ —ipsiwd — —ik\\ (cf. Fig. [3). 
Thus, closing the contour in the lower kz plane picks up 
two residues, one at — and one at —ipsin-d. The result 
is 
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e 

8^ 



e 

4^ 



dp 



dt — 

-1 P- 



2ipot 



7|| -Pll -7oip< + 



TM 



it) e 



2ap 



dp / dt 



p2 + 4p2t2 



i2 



7|| 



t^ 

- 270J50 (2 



7|| 



t^ 



t^) + 



7|| • Pll 



7|| ■ P|l - 27oipi + 273P 
Pll - 2-foipt + i?^M(0) e-2-P^/T^ 

l-i2 



i2 



-273P2 



TM 



it) e 



-2ap 



47oPor + 273pJi?TA/(0)e 



.(43) 



where we have again renamed cos -d = t and abbreviated 



- ^{n^ - l)t^ + 1 



As before, we are interested in an asymptotic result for E^*^ for large values of poa. To be able to do asymptotic 
analysis one needs to separate the terms with different arguments in the exponential. However, doing this simple- 
mindedly leads to two divergent integrals because their integrands each behave as 0{t^^) for t ^ 0. That is why we 
add and subtract the same term and subdivide the integral as follows. 



(44) 



with 



-'C 



Zjt-I — 



e 

4^ 
47r2 

o2 



d^ / dt 



(27 



7|l - PI 
7n - Pi 



- [270P0 {t^ - 2) - 273P.] R^^Mit) e-^P"'^*^ 



d^ / dt 



47r2 



47oPoi?TA/(0) 



R 
d^ 





TA/(0) 
1 







^2 + U2 
d^ dt- 



IR"^ 



■TM 



it) 



TM 



(0)] 



— 4poatvT— t^C 







1 



where we have again rescaled p = 2pQt^. 

We do the asymptotic analysis of these integrals one by one, starting with S^*^. In order to get an asymptotic 
expansion for large pga, one would try to integrate by parts. However, The integrand of Sj^^ behaves as 0(1) 
for t — > 0, and thus the factor 1/t that one gets through integrating by parts in the ^ integral would destroy the 
convergence at i = 0. Adapting the general method of obtaining an asymptotic expansion of such two-dimensional 
integrals ^-dd and subtract the problematic point at i = and write 

r-l 



dt 



1 



(47oPo + 273P,)i?|^M(0) e-4P""*« 



dt 



- { [270P0 (i' - 2) - 273P.] Ri^Mit) + (470P0 + 273P.) i?^A./(0)} e-^f'"^'^ 



The first of the integrals is easy to calculate; the t integral is immediate, and the remaining integral over ^ is a well 
known combination of sine and cosine integrals [l3|- The integrand of the second integral now behaves as 0{t^) for 
t ^ 0, and we can thus integrate by parts twice without getting convergence problems at f = 0. In this way we find 
to order l/(poa)^) 



^A 



Stt^Poq 



(27oPo+73P.)-Rta/(0) 



1 



4poa 



327r2p2a2 



{loPo ^ ^ [it' - 2) RffMit) + 2 R^iiO)] 



'^TMit) 



^TA/(0)] 



(45) 



r 



for which we have also made use of the known asymp- 
totics of the sine and cosine integrals [T^. 



The asymptotics of can be calculated similarly 

by adding and subtracting the next term in the Taylor 
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expansion of R^^,i{t), i.e. by replacing 



^2 



The first term can then be integrated exactly, and the 
rest behaves as 0{t^) for t — > and can thus integrated 
by parts with respect to ^ twice. The result to order 
l/(poa)^) is 



— hlaPz 



-711 - Pi 



1 ( 1 



^R 



L II 

■TM 



(0) 77 



1 



1 



dt 



- ^R 



TM 



'(0) 

. 

Next we turn our attention to the asymptotic evaluation of Sp^^. We cannot separate the two summands in the 
integrand because otherwise the t integral does not converge. Thus, to manipulate just one part, we must set the 
lower limit of the t integral to some small positive s and take the limit e — > only once we have combined all parts 
again. We write 



at 7: e 



r 



1 

1/72 



In the integral that runs from l/\/2 to 1 we make a change of variable from i to s = — t'^. Then renaming s into 
t again and ignoring terms that vanish in the limit e — > 0, we find 



at ( 



dt 



+ 



(1-^2)3/2 



Now we substitute v — t\/l — t^ and obtain 



at ^ e 



1/2 



dv- 



l-2z;2 



-4^00-0^ 



^2^1 - Av"^ 



Renaming the integration variable v into t again, we can use this identity to write S^*^ as 



e2 / V f /"^/^ 1 / 



l-2t2 
VI - Af^ 



-4po at^ 



1/2 



dt-e-^P""^^ 
t^ 



Integrating by parts in the integral twice, we obtain to 
order l/(poa)^) 



e 

1^ 



+7|i •PIi)^tm(O) 



2^2 



(47) 



obtain 



1/2 

r7oPoi?TM(0) \ di \ dv 



-^IoPoRtmC^) 



/•1/2 
d( / dv 
Jo 



1 



$2 + 1 Vl - 4i;2 



.1/2 
dU dv- 
Jo ? 



X e 



1 



-1 e 



-4poa< 



For the asymptotic expansion of ^J/^ we apply much 
the same tricks as above: we split the t integration at 
l/-\/2 and in the integral from l/-\/2 to 1 we substitute 
s = Then we rename s back to t, combine 

the two integrals again, and substitute v = t\/l — t^ to 



The first of those integrals can be solved in terms of 
known special functions |1J|, and in the second integral 
we can integrate by parts twice with respect to ^. Thus 
we obtain for S'^*'^ to order l/(poa)^) 



(48) 
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According to (|44|l we combine the results H45H48|I to the TM modes. It turns out to be one order larger than 
obtain the leading-order self-energy contribution from that from the TE modes. 



327rpoa 



-n? (it? — l) 



7|| • Pll 



(n2 + 1)^ 



273^2 



- 



(n2 + ly 



O (1/poa) 



Thus, to leading order the radiative self energy is due to just the TM modes. 



327rpoa 



7|| • Pll 



[n' - 1) 
(n2 + 1)^ 



273Pz 



2n^ - _ 1 
(n2 + if 



(49) 



The next-to- leading term is easily determined from the results H42I) and H45H48|I . While all the t integrals are 

elementary, it is convenient to use formula manipulation software like Maple to evaluate and combine them. The final 
result for the next-to-leading order contribution to the radiative self-energy is 



327r2(poa)2 



270P0 



1 



In ( n + \/ n? — \ 



n? + 4n2 + n^2 



2^4 



-7n -P 



n2(2n4 + 3„2 _ 3^+ 1) 



(n2 -f 1) (n -I- 1) (n2 + l) 



3/2 



arctanh- 



n- 1 



2^4 



3(n+l)(7i2 + l)2 (n2 - 1) (n2 + 1)^/2 



arctanh- 



3(n-h l)(n2 + 1)2 



(n2 - 1) (n2 + 1)5/2 



rt — 1 
" %/n2 + 1 



arctanh- 



n — 1 
+ 1 



D. Self energy in the limit n ^ 00 

In the limit of perfect reflectivity n — > cxd the calcu- 
lation of the self energy simplifies considerably. All re- 
flection coefficients go to either -1-1 or -1 (cf. Eq. IjAlOjl '). 
and the photon propagator takes on a much simpler form 
with the fcs integration running straight along the real 
axis (cf. Eq. dJ). The calculation of the self energy 
can then proceed in exactly the same way as explained 
in section IIV Al above. It starts to differ only with the 
asymptotic analysis of Eqs. (|1T|| and For TE the 

asymptotic analysis of 1)41(1 relied on the fact that R^^ (t) 
behaves as 0{t^) for t 0, but in the limit n ^ 00 we 
have = — 1 for all i, which leads to a very different 
asymptotic behaviour of the self energy E-^-^ . To leading 
order we find 



yTE , 
perf ' 



327rpoa 



7|| 



In the case of TM, something similar happens. The inte- 
grals Sp*^, Sl^*^, and to leading order in l/poa also S^^*^ 
give the same with the limit n — > 00 taken first as they do 
for finite n and with the limit n — > 00 taken in the end re- 
sult of the asymptotic calculation. However, the integral 
S^*^ does not even appear if the limit n — > c» is taken 
straightaway. For finite n its asymptotically leading term 
depends on the second derivative i?y^j"(0), which is of 
course zero if the limit n — > cx) has been taken first and 



Rx]^ is a constant. If one takes n 



first then to 



leading order only Sj*^ and S-^*^ contribute, and one 
obtains 



yTM , 
perf ' 



167rpoa 



In total the radiative part of the self energy with the 
perfect-reflector limit taken first is 



Sper/(P, -2a) 



327rp( 



7ll • Pll - 273PJ 



which clearly differs from the limit n ^ 00 of Eq. 149|l . 
So, mathematically not surprisingly, we find that the re- 
sult for the self energy differs depending on whether we 
perform the calculation for finite n and subsequently take 
the limit n —f 00, ov whether we take the limit of per- 
fect reflectivity n — s- cx) flrst and evaluate the integrals 
then. The cause of this is simply the fact that the limits 
n — > 00 and t — > in the Fresnel coefficients R^ (t) do not 
commute. The point t — Q corresponds to cos-d — 0, i.e. 
to /co — 0. Thus, physically speaking, the limit of perfect 
reflectivity is not interchangeable with the limit — > 0; 
in other words, long-wavelengths excitations get treated 
very differently in and away from the limit of perfect re- 
flectivity, n — > 00. Indications of this can also be seen in 
the non-relativistic calculation of the energy level shift 
'7]. 
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The important lesson to be learnt from this observa- 
tion is that models that assume perfect reflectivity from 
the outset are bound to give the wrong answer if long- 
wavelengths excitations play any role in the system under 
investigation. Luckily, most of cavity QED is concerned 
with atoms and other bound systems which have an in- 
herent low-frequency cut-off (e.g. the lowest transition 
frequency of an atom in state \i) to dipole- allowed 
states \j)). However, for unbound or partially bound sys- 
tems a perfect-reflector model is principally inadequate 
for describing any physically realizable system, no matter 
how good the reflectivity of the boundaries may be ■ 

E. Electrostatic contribution 

The evaluation of the electrostatic shift H29() is straight- 
forward. The state in which the expectation value is be- 



ing taken is a wave packet that is localized at around 
= (0, 0, —a) and that has an average momentum p. 
Here the localization parallel to the surface of the di- 
electric can of course be arbitrarily loose, as the system 
is translation invariant parallel to the surface and hence 
the energy shift does not depend on the transverse loca- 
tion of the particle. We choose to represent the localized 
state by a Gaussian wave packet, 



Taking the expectation value in Eq. (|29|) in this state and 
using the canonical mode expansion H18|l for the spinor 
operators, we obtain 



xe-^'^-^^^ J2 u^'\q') 7° u^^\p') u'^^\p') 7° ^('^(q) G<i.(x,x') . 
i=i,2 



Since the integrand as function of x and x' is peaked at 
the location x^ of the wave packet, we can approximate 
the Green's function G$(x,x') by G'$(xa,Xa). Then the 
X and x' integrations are easy to carry out and give 5 
functions. The sum over polarizations j can also be done 
because the Dirac eigenspinors satisfy 



J2 u^^\p')u'-^\p') 
i=i,2 



1 
2m 



m) 



Thus the above expression simplifies to 

^2 



5E, 



Coulomb 



G<i.(x„,xJ / d-'p' — 



m (p -p) 



2^3/2^3 - J . 2p';^ 

X (p')(7°Po + 7 • P' + my/'Hp') ■ 

The integrand of this expression peaks at p' = p, so that 
we can approximate p' by p everywhere except in the 
exponential and carry out the integration. Taking into 
account that p is on shell, we then find 



2 

~ ^ G<I>(Xa,Xa) — ?l(')(p)7°u(')(p) 

2 Po 



It remains the evaluation of the Green's function at 
the location of the particle. Since we are looking for the 
energy shift relative to a particle in free space, we use the 
difference of the Coulomb Green's function (|16|l and the 



free-space Coulomb Green's function 



G^{x, x') = G^{x, x') - j e'(--') k J_ _ 

As Xa is outside the dielectric, we have 

G4,(Xa,X,) =--^^— / / 2 



n'^ + 1 47r2 Jo 



For the integration we close the contour in the lower 
half-plane. The integration over fcy is then trivial. The 
result is 



G$(Xa, Xa) — — 



-1 1 
v? + 1 Stto 



Thus the electrostatic energy shift is 
'n? — 1 m 



5E, 



Coulomb 



uW(p)7%(*)(p) . (50) 



For a particle at rest this agrees with the classical energy 
shift of a point particle in front of a dielectric half-space 

If one wishes, one can express this shift as being due 
to a Coulomb self-energy function. If the shift is given 
by Eq. (|28|l then we can write 



^Coulomb 

(p, -2a) 



n'^ + 1 Wna 



7" • (51) 



15 



SUMMARY AND DISCUSSION OF THE 
RESULTS 



the total self-energy to leading order in l/{poa) 



Combining our results for the radiative self-energy H49|l 
and for the electrostatic self-energy (|51|) . we obtain for 



■Jtotal 



(p, -2a) 



32npoa 



7|| • PI 



n2 - 1) 
(n2 + if 



2j3Pz 



- 



1)^ 



27%- 



(52) 



The total energy shift is easily determined from Ea. (|28|) . Since spin-up and spin-down states are degenerate without 
the perturbation, the right-hand side of Ea. H28|) is actually a matrix with the energy shifts as eigenvalues. In general, 
the self-energy operator (|52|l is not diagonal in the spin states ^''•'(p), as given in Eq. H19() . if spin up and down states 
are defined along z, i.e. for 



a(i) 



(2) 



We find 



'^-(i)v 
hi 



(1) = -Il(2)i]. 
E 



total" 



(2) 



327ra£;2 



and 



m 



U^^JStotalU 



(1) 



Thus the energy shift is 
SE 



total" 



(2) 



' n'^{n? — 1) 



n^jn' - 1) 
(n2 + 1) 



1 E 



1 E + m 



'S2TraE^ 



2n^ — — 1 , n^fn^ 



1 



(n2 + 1)2 



(n2 1) 



n2 - 1 



-E' 



± 



n2(n2 — 1) 
(n2 + 1)2 



n2 - 1 

) 

n2 + 1 _B + m 

I 



(PxPz)'^ + {PyVzY 



-E' 



{{Px + iPy)Pz 



(53) 



(54) 



(55) 



For wave packets that are either stationary or whose mo- 
tion preserves the symmetry of the problem, the non- 
diagonal elements H54|) are zero and the energy shift is 
given simply by (|53|l . 

In the limit of perfect reflectivity n ^ oo the calcula- 
tion yields a result for the total self-energy that differs 
from the limit n —f oo of Eq. H52I) , 



^perf 
■^total 



(p, -2a) 



32TTpoa 



-7|| • PI 



-t- 2j3pz + 27% 



(56) 

Accordingly, the energy shift differs from the limit n 
oo of Eq. (|55|l . We have discussed the reasons for this 
discrepancy in Section llV Dl and in Ref. j). 

Quite apart from calculating the energy shift of an elec- 
tron in front of an imperfectly reflecting half-space, we 
have established the major building blocks for QED in 
the presence of a dielectric half-space. Two alternative 
formulations for the Feynman propagator of the electro- 
magnetic field are given in Eqs. H1U|I and Hll|) . The loss of 



translation invariance perpendicular to the surface of the 
dielectric half-space is an essential complication in loop 
calculations, but we have demonstrated how to tackle 
this at one-loop level. 
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APPENDIX A: POLARIZATION VECTORS AND 
NORMAL MODES 



In generalized Coulomb gauge the direction of the elec- 
tromagnetic field can be described by the following choice 
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of polarization vectors: 



For the right-incident modes the vector potential reads 



= -Gu G 



= (-A||)-i/2(0,-z9„ia,,0) , 



G - -^t^ G = (-A) ^/^(O, -id^, -idy, ~id^) , 
e'^ = e^c = (l,0,0,0). (Al) 

Here A ~ d'^+dy+d'^ is the Laplacian in three dimensions 
and A|| = d^+dy is the one in two dimensions parallel to 
the surface of the dielectric. The physical polarizations 
are the transverse electric and the transverse mag- 
netic e^jyj, which have vanishing electric or magnetic, 
respectively, field components perpendicular to the sur- 
face. The unphysical polarizations are the longitudinal 
and the timelike e'^. 

Constructing the normal modes is straightforward if 
one proceeds from a plane incident wave which gets re- 
flected and transmitted at the surface. The incident wave 
is what it would be in free space (for left-incident modes) 
or in a homogeneous dielectric (for right-incident modes), 
and the transmitted and reflected components can be de- 
rived from the continuity conditions (0). 

For the vector potential of the left-incident mode one 
finds 



a=TE,TM,G,C 



(27r)3/2(2cj)i/2 

x<{e^exp(-ia;°A:°)a^(A:) 

-z) (exp(ik • r) + R^{k) exp(ik'^ • r)) 



-F6i(z)-T^(fc)exp(ik''-r) +h.c. ^, (A2) 
n J J 

where a^{k) is the photon annihilation operator of the 
mode and the reflection and transmission coefficients are 



Rte{^) 



h — h'^ 

k^ + k<i 



2nk, 



h 4- hd ^ 



n'^kz + kz 



n 



21, _|_ l.d 



TSik) 



- n^ki 



kz - n^kj 
kz + rfik'^^ ' 

Th{k) = 



2nk, 



(A3) 



The wave vectors are 

k = (k|| , kz) , k"^ = (k|| , -kz) , = (k|| , kf) , 

ki = sgn{kz){n'^kl + {n^ - l)fcf )i/2 , kj = kl + kl , 
ko = u;, iv^^kl + k^ = \{ki^ + kh. (A4) 

II 7-j II 



E 



d^k e{-ki) 

(27r)3/2(2cj)i/2 



(T=TE,TM,G,C 

x|e^exp(-ia;°A:°)a^(A;) 
1 



)- (exp(zk'' • r) + R^{k) exp(ik'^''* • r)) 



+9{~z)T^{k) exp(ik • r) -I- h.c 



with the reflection and transmission coefficients 



(A5) 



R k"^ — kz 

Rrp^[k) = , ^ , , , Trp'^{k) 



2fc2 



R. 



■TM 



z+kd 

(k) = ir-^ , T^Ak) 



n{k^ + kz) 
2k'i 



k'} + n^kz 



kf + 11? kz 



Rc{k) = R^ik) 



2i,d 



n'^k'^ + kz 



T^ik) 



2i,d 1 i. ' ^G (fc) - „2;jd .r 1. ■ (^^) 



n^k 



n 



Because these modes are right-incident the z component 
of the incident wave vector is negative kf < 0, and the re- 
fiected wave has the wave vector k'*''' = (k", —kz). Note 
that the integration over k'' includes imaginary values 
of kz , which correspond to modes that come from inside 
dielectric, suffer total internal reflection at the interface, 
and are evanescent on the vacuum side. One has 



kz = sgn(fcf)-./fcf - (n2 - 1)A:2 

n V II 

for fcf 2 _ („2 _ i)fc2 > 



1 



= -J-y'-fc^2 + („2_l);.2 

for kf - (n^ - l)fc2 < 



We have chosen the branch of the square root such that 
the evanescent modes are truly evanescent, i.e. exponen- 
tially falling away from the interface on the vacuum side. 
This also ensures that these modes are genuinely totally 
reflected, i.e. that the relation R^R^ = 1 is fulfilled. 

For the physical polarizations TE and TM the above 
modes are well-known jl^ . It is easy to see that they are 
mutually orthogonal, but the proof of their completeness 
is surprisingly tricky [l5l |. The following relations are 
useful for showing the completeness of the modes and 
for simplifying our expressions in section FlI Al They are 
valid for all polarizations so that we drop the index a. 
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For real k, we have 



[TLT^*]{k,,kf) + [R^R^*]{-k,, -ki) = 1 



i?^*(-fcf) = i?«(Af) 



n^k 

and for imaginary k. 



(A7) 



= 1 



i?«*(-fcf) = i?«(fcf) 

[R^*T^]{-kf)^T^{kf) 

[R^T^*]{-ki)^T^*{ki) 



(A8) 



Applying the polarization vectors HA1(I on a plane wave 
e , as m (|X2)l and (|X5)l . one can express them in terms 
of the wave vector k. However, it is important to realize 



that the incident, transmitted, and reflected components 
all have different wave vectors and thus, according to 
(lAip . have polarization vectors that point in different 
directions. All four polarization vectors form a complete 
and orthogonal system, and the TE and TM polarizations 
are complete in the subspace of physical states 



cr=TE,TM,C,G 

J2 = g^"^ - ry^r;^ + k^k'' . (A9) 

a=TE.TM 



Here k'^ ~ (k^ ~~ {k'r\)'rf') j ^ [kr])'^ — fc^ is the unit vector 
along the space-like part of fc'', and r/^ = (1, 0, 0, 0). 

Finally, we would like to consider the limit n — > oo 
which is commonly thought of as corresponding to a half 
space bounded by a perfectly reflecting wall. Indeed, 
in this limit the reflection and transmission coefficients 
become 

R^{k) = R^{k) ^ -1, i?^s(fc)^-l, R^j,jik)^l, 



-R^ik) ^ , -T^^ik) ^ , rj^(fc) 



0, 



(AlO) 



so that only the left-incident mode l|A2|l survives and gets 
perfectly reflected at z = 0. 
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